Analysis of resonant inelastic x-ray scattering in La 2 Cu04 
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We provide a semiquantitative explanation of a recent experiment on the resonant inelastic x-ray 
scattering (RIXS) in insulating cuprate La2Cu04 (Y.J. Kim et al., Phys. Rev. Lett. 89, 177003 
(2002).). We show theoretically that there are three characteristic peaks in RIXS spectra, two of 
which are attributed to the charge transfer excitation and are reasonably assigned to those observed 
experimentally. The lowest energy peak has a relatively large dispersion (~ 0.8 eV) and is suppressed 
near the zone corner (7r,7r), in agreement with the experiment. We stress that electron correlation 
is an essential factor for explaining the overall energy-momentum dependence of the RIXS spectra 
consistently. 

PACS numbers: 78.70.Ck, 74.72.Dn, 78.20.Bh 



I. INTRODUCTION 

Recently it has been revealed that resonant inelastic 
x-ray scattering (RIXS) is a powerful tool for elucidat- 
ing the electronic properties of solids 1 . RIXS is a unique 
technique of detecting the relatively high-energy (of the 
order of 1 eV) charge excitations in solids. It is naturally 
expected that excitation spectra in such a high-energy 
regime reflect not only detailed electronic structures over 
a wide energy range but also electron correlations origi- 
nating from strong electron-electron Coulomb repulsion. 
However, unfortunately, there are still only a few theoret- 
ical works on analyzing the effects of electron correlations 
on such high-energy excitation properties. 

RIXS measurement in the hard x-ray regime has been 
applied to search for charge excitation modes and to de- 
termine the momentum dependence of charge excitation 
in several transition metal oxides, with the incident pho- 
ton energy tuned to the absorption K edgeSi^^. A 
number of insulating cuprates have been investigated up 
to now, under the expectation that such a kind of exci- 
tations might be related to the unconventional high-T c 
superconductivity. For example, Kim et al. recently re- 
ported a detailed study of the momentum dependence of 
the charge transfer (CT) excitation spectra in La2Cu04 
(Fig. P; 5 . They observed two characteristic peaks in the 
spectra at around ui ~ 2.2 eV and ~ 3.9 eV at the 
zone center, by setting the incident photon energy to 
the Cu absorption K edge. The low-energy peak at 2.2 
eV shifts to u> ~ 3 eV around (tt, 0), and becomes very 
small around the zone corner (tt, 7r). They considered 
that the two peaks with such a highly dispersive behav- 
ior are attributable to some kind of bound excitons which 
are formed as a result of Coulomb interaction, since such 
a bound exciton could move without disturbing the an- 
tifcrromagnetically ordered background. 

In this article, we provide an explanation of the RIXS 
spectra for the typical antiferromagnetic charge transfer 
insulator La2Cu04 at a microscopic scale, by calculating 




FIG. 1: RIXS spectra obtained by Kim et al. are shown for 
various momentum transfers as a function of energy transfer. 
(Reproduced from Ref. 0.) 



the transition probability in the RIXS process. In or- 
der to investigate electron correlation effects in the RIXS 
process, we use the perturbative method developed by 
Platzman and Isaacsi The scattering probability is cal- 
culated in the manner adopted by Nozieres and Abra- 
hams to discuss Fermi edge singularity in metals^. We 
use the Hartree-Fock (HF) theory to describe the an- 
tiferromagnetic ground state, and take account of the 
electron correlations in the scattering process within the 
random phase approximation (RPA). The perturbative 
method for band models compares well with other theo- 
retical techniques based on a local atomic picture, such 
as exact diagonalization for a small cluster—. We obtain 
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a three-peak structure in the spectra at around 2, 5 and 
9 eV, by taking account of the 02p band explicitly. The 
spectrum obtained as a function of energy loss shows dis- 
persive behavior depending on the momentum transfer. 
In particular, the lowest energy spectral weight around 
2 eV shows a relatively large dispersion (~ 0.8 eV), and 
becomes small around (7r,7r), in agreement with the ex- 
periment. The two peaks at around 2 and 5 eV corre- 
spond to charge transfer excitations and are assigned to 
those observed in the experiment by Kim et al£. 

II. MODEL AND THEORY 

We introduce the model Hamiltonian for the system in 
the form H = H dp + Hi s _ 3 d + H ls + H 4p + H x . H dp , H ls 



where dk a and p-kio {d k(T anc ^ PLea) are annihilation 
(creation) operators for Cu3d x 2_ y 2 and 02pi electrons 
(£ — x,y) with momentum k and spin a, respec- 
tively. We adopt the dispersion relations £d pe (k) = 
2it dp sin !f, ( ptpe (k) = e p = 0, and ( Pj;Py (k) = C PyPx (k) = 
it pp sin 4f sin , where t dp = 1.3 eV and t pp — 0.65 
eV, which were estimated from LDA calculations^. The 
Coulomb energy and the charge transfer energy in HF 
theory are Udd = 11 eV and e^ F = —0.7 eV, re- 
spectively. Within HF theory, we obtain the antifer- 
romagnetic ground state with staggered spin moment 
m = 0.46/iB for the present parameter set. H\ s —3d gives 
the interaction between Culs holes and Cu3g? electrons: 

H ls - Sd = — ^ 4'+qa4-q CT ' S k -'eZkV, ( 2 ) 
kk'qcro"' 

where V is the core hole potential and is taken to be 
15 eV, and Sk<r (s fci is the annihilation (creation) oper- 
ator for Culs electrons with momentum k and spin a. 
H x describes the transitions between the Culs and Cu4p 
states, which involve photon absorption and inverse emis- 
sion processes. H x is of the form 

H x = ^( W ( q;k K T +qa S ka + h.C.), (3) 
kqtr 

where p 1 * is the creation operator of the Cu4p electron. 
Since the Culs wave function is highly localized in the 
coordinate space, the wave functions of the Cu4p band 
states only in the region close to the origin contribute to 



and 7?4 P describe the kinetics of Cu3d :c 2_j / 2 (hybridized 
with 02p orbitals), Culs and Cu4p electrons, respec- 
tively. We take a completely flat band for the Culs elec- 
tron, which is well localized in the solid. In addition, 
for the Cu4p electrons, we take a simple cosine-shaped 
band with the energy minimum at the zone center. In 
any case, no detailed structure of the Cu4p band is re- 
quired to determine the shape of RIXS spectra, since the 
factor including the dispersion function of the Cu4p band 
is integrated up with respect to the momentum and does 
not affect the momentum- energy dependence of the spec- 
tra, as we see later. For H dp , we use a two-dimensional 
dp-modeh 



(1) 




FIG. 2: (a) Diagram for the transition probability in the 
RIXS process. Green's functions for Culs, 4p and 3d elec- 
trons are assigned to the solid lines with 'Is', '4p' and '3d', 
respectively. The wavy and broken lines represent the pho- 
ton propagator and core hole potential V, respectively. The 
shaded triangle is the effective scattering vertex renormal- 
ized by Coulomb interaction Udd- (b) Vertex renormalization 
in RPA. The empty square represents antisymmetrized bare 
Coulomb interaction between Cu3d electrons. 

the matrix elements. This strong local nature results in 
the momentum independence of the matrix elements. In 
the present article, we ignore the momentum dependence 
of the matrix elements u>(q; k), i.e., w(q; k) = w. 

In the RIXS process in La2Cu04, the following scat- 
tering channel is expected to predominantly contribute 
to the main spectral weight: the incident x-ray with the 



Hd v = e ddl a d ka + C PePe , 0^pt l(7 Pke'a + J2(£d Pe (k)dl a p kea + h.c.) 
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ls-Ap absorption energy excites the Culs electron to the 
Cu4p state. In the intermediate state, the created Culs 
core hole plays the role of a localized strong scattering 
body for the Cu3d electrons near the chemical poten- 
tial. After the core hole scatters the Cu3d electron sys- 
tem from the ground state to some excited states (the 
most probable scattering process is the creation of only 
one electron- hole pair), the core hole is annihilated to- 
gether with the Cu4p electron and a photon is emitted 
in the final state. Compared with the incident absorbed 
photon, the emitted photon loses momentum and energy 



in amounts equal to those transferred to the Cu3d-02p 
electrons. In the final state, after the Culs core hole dis- 
appears, an electron-hole pair remains on the Cu3d-02p 
band. In the present study, we consider the lowest order 
in the core hole potential V, i.e., simple Born scatter- 
ing of Cu3d electrons by the core hole. The transition 
probability for the total process is diagrammatically ex- 
pressed in Fig. |2I a) ■ The diagram in Fig. ffia) is con- 
verted straightforwardly to an analytic expression in the 
same way as described in Ref. H: 



J 



W( qi , q f ) = (2ir) 3 N\w\ 4 £ d(E,(k) + uj - E f (k + q)K(k)(l - n r (k + q)) 



CWk)A CTCT ^, q)Ul, tjl (k + q) £ 



V 



ki 



(uji + ei s + iT ls - e 4p (ki))(w/ + e ls + iT ls - e 4p (ki)) 



(4) 



where characterizes the diagonalized band, Ej(k) 
is the energy dispersion of the band j, qi — (u>i,qi) 
[qj = (wj,qj)] is the energy-momentum of the ini- 
tially absorbed [finally emitted] photon, and the energy 
and momentum transfers are given by lo = u>i — u>f 
and q = q» — q/- eis and 64 p (ki) are the kinetic 
energies of the Culs and CuAp electrons, respectively. 
rij(k) is the electron occupation number at momentum 
k in energy band j. [7j i( jo-(k) is the (j, da)- element 
of the diagonalization matrix of the HF Hamiltonian. 
We set the incident photon energy near the absorp- 
tion edge, i.e., u>i « 6^(0) — e ls , and assume that 
the decay rate of the Culs hole is Ti s — 0.8 eV. In 
the present study, we simply use the vertex function 
h-atji (q) within RPA. The vertex renormalization within 
RPA is diagrammatically represented in Fig. |2Jb). It 
is interesting to note that the factor nj(k)\Uj y da(k-)\ 2 
[(1 - 7v(k + q))|t/j CT/j .,(k + q)| 2 ] in Eq. Q represents 
the partial occupation number of the Cu3d electron [hole] 
with momentum k [k + q] and spin a [a 1 ] in the band 
j [j']. According to our investigations, the product 
r^(k)(l-rcj/ (k + q))|C/ 3 - da (k)[/] CT , -,(k + q)| 2 is the dom- 
inant factor in the scattering amplitude. Therefore, the 
RIXS spectrum in La2Cu04 reflects the Cu3<i partial 
electron and hole occupation numbers, that is, the Cu3d 
partial density of states. 



III. RESULTS 

The numerical results of the RIXS intensity calculated 
using Eq. (0J are displayed in Fig. |3] We see the three 
characteristic peaks at around u> ~ 2, 5 and 9 eV. The 
lowest energy peak at 2 eV shows relatively large disper- 
sive behavior. The peak shifts by approximately 0.8 eV 



at q « (it, 0) and becomes small at q « (k,tt). This 
characteristic momentum dependence is consistent with 
the results in Fig.^ The peak at 5 eV is assigned to the 
peak observed experimentally at 3.9 eV. The crude tight- 
binding fitting would be responsible for the 1 eV devia- 
tion of the peak position, and would still be insufficient to 
reproduce the electronic structure in such a high-energy 
region. The two peaks at 2 and 5 eV originate from 
CT excitation: Cu3d electrons hybridized with the 02p 
band are excited to the upper Cu3d band. This two-peak 
structure is a result of the structure of the Cu3c? partial 
density of states mixed in the broad 02p band. On the 
other hand, the small peak at around 9 eV corresponds to 
the transition from the Cu3cZ lower Hubbard band to the 
CuSd upper Hubbard band. The transition probability 
for this 9 eV process is relatively low. This is because the 
majority spin state at a local Cu site is almost fully occu- 
pied by one electron, and the hole occupation number in 
that state is almost zero. The actual localized moment 
will be almost 1(J-b, since the electron occupation num- 
ber per one Cu site is almost one. As we have mentioned 
above, the intensity is closely related to the product of 
the Cu3<i electron occupation number and the Cu3<i hole 
occupation number. In the local picture, this product is 
expressed as n,j j j (T (l — rid.io-) (jid,i<j is Cu3c£ electron num- 
ber at Cu site i with spin a) , and becomes almost zero in 
the limit of m = 1/ib- In the present HF calculation, we 
may have underestimated the electron occupation num- 
ber in the majority spin state at a Cu site. In other 
words the magnetic moment m = 0.46/iB obtained in 
the HF theory might be smaller than the actual localized 
moment. Therefore, we should consider that the present 
method still overestimates the transition probability for 
this 9 eV process. Actually, the intensity for this process 
would be too small to be observed. Therefore this 9 eV 
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FIG. 3: The upper figures show the Cu-0 plaquette and 
the contour in the first Brillouin zone (BZ). The lower six 
panels show the RIXS spectra as a function of the energy 
loss lu = uii — ujf on the above contour in BZ. The thick solid 
(thin broken) curves are results obtained with (without) RPA 
vertex renormalization. 



weight should not be identified with the 7.2 eV peak ob- 
served in the experiment (peak c in Ref. 0), whose origin 
should be discussed on the basis of other scenarios. 

Here we inspect the electronic structure and the den- 
sity of states (DOS) of the model in the antiferromagnetic 
ground state. The Cu3d partial density of states (PDOS) 
and the total density of sates are displayed in Fig. Q] The 
2 and 5 eV RIXS peaks are assigned to the electronic ex- 
citation processes from the Cu3d PDOS weights around 
lu w — 1 and —4 (eV) to the upper Cu3g? band around 
to w 2 (eV), respectively. These Cu3g? PDOS weights 
around to sa — 1 and —4 (eV) originate from the strong 
hybridization of the Cu3g? and 02p orbitals. Therefore 
the two peak structures are attributed to Cu3d-02p CT 
excitations. 

It is of interest to discuss the effect of the correla- 
tions on the shape of the spectra. In Fig. [21 we com- 
pare the RPA-corrected case with the uncorrected case 
(i.e., zero order in Udd- J^aa'(q) 5 aa r in Eq. Q). We 
find that the lowest order uncorrected calculation repro- 
duces roughly the overall structure of the RIXS spec- 
tra. This is because the spectral weight is determined 
mainly by the partial occupation number of Cu3g? elec- 
trons in each band, and the detailed scattering process 
is not necessary in the first step of analysis. However, 
the low energy peak structure at around 2 eV is modi- 




-4 -2 
co (eV) 

FIG. 4: The total density of states and the Cu3d x 2 _ y i partial 
density of states are shown by the solid and broken lines, 
respectively. The inset shows the band dispersion of the model 
in HF theory. 
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FIG. 5: Left panel shows details of the spectra along the line 
(0, 0)-(-7r, 0) between lu = 2-3 eV. The data are offset vertically, 
and the thick horizontal bars drawn on the left vertical axis 
denote the baseline for each spectrum. The arrows denote the 
positions of the main peak. Right panel shows the position of 
the main peak as a function of q. 



fied somewhat by the RPA corrections. The correlation 
effects between Cu3d electrons cause an anisotropy and 
modify the simple s-wave scattering amplitude, and re- 
sult in the enhancement and suppression of the spectral 
weight at q = (0, 0) and q = (it, 7r), respectively. 

We focus our attention on the dispersive behavior of 
the low-energy peak at 2 eV. We display the detailed 
behaviors of the 2 eV peak in Fig. EI along the line (0, 0)- 
("7r, 0). We find that there is a fine structure in the 2 
eV peak: the peak splits into two peaks. Such a fine 
structure is not observed in the experiment—, maybe due 
to the low resolution. The shift of the peak between the 
points (0,0) and (n, 0) is about 0.8 eV. This is slightly 
smaller than the experimental result. We consider that 
this deviation originates from our underestimation of the 
02p bandwidth. 
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IV. DISCUSSION 

We have used the LDA parameters for the model. Al- 
though there are still much room for discussing the valid- 
ity of the application of the LDA parameters to excitation 
processes, the parameters would provide a good starting 
point. To be more quantitative, the renormalization ef- 
fects due to correlations should be taken into account in 
the future. 

We briefly note some effects of other factors which 
are not included in the above calculations: effects of 
on-site Coulomb interaction U pp between 02p electrons, 
inter-site Coulomb interaction Ua P between 02p and 
Cu3d electrons, higher orders in the core hole poten- 
tial V, self-energy corrections and Coulomb interaction 
^3d-4p between Cu3g? and 4p electrons, multi-pair cre- 
ation, phonons. We expect that these effects are not 
crucial in explaining the observed spectral properties. 

We employed multiband RPA including U pp , but found 
that U pp only reinforces the effect of the Coulomb in- 
teraction Udd- We cannot exclude the possibility that 
Udp modifies the analytic properties of the spectral func- 
tions, i.e., novel poles of the spectral function could ap- 
pear on the complex to plane near the real axis. In such 
a case, some peak structures related to a bound exciton 
mode would be added to the spectra. However, note that 
there is no necessity for the additional peak to be located 
around the observed positions at 2 and 4 eV. Such a novel 
structure would not be consistent with the experiment. 

We could investigate the effects of higher orders of V by 
summing up the perturbation terms to the infinite order, 
similarly to T-matrix approximation for single impurity 
problems. We speculate from our preliminary work that 
the higher orders in V only negligibly modify the spectral 
shape, although they modify the absolute intensity of the 
spectra. 

Concerning the self-energy corrections, detailed inves- 
tigations are a future work. In x-ray absorption process, 
U^d-Ap could play an important role for "Extrinsic losses" 
(e.g., Ref. Illf) through self-energy corrections, leaving 



electron-hole pairs on the 3d band in the final state. But, 
U?,d-ip would not be so essential in the present RIXS 
process, since the Ap electron is annihilated rapidly with 
Culs hole in the final state of the RIXS. 

Multi-pair creation processes may become more impor- 
tant for metals, leading to the Fermi edge singularity in 
the spectra. However multi-pair creation processes are 
not so essential in insulators, since they require much ex- 
citation energy (the energy per one pair is given by the 
insulating gap, about 2eV in the present case). 

Regarding phonon effect, phonons can be excited in 
order to screen the core hole potential in the intermediate 
state. However, in the present RIXS, the phonon process 
cannot be the dominant process to screen the core hole 
potential, since the annihilation of the core hole is much 
faster than the screening process by phonons. 

V. CONCLUSION 

In conclusion, we note the following points, (i) We 
have not employed the physical picture that a kind of 
bound exciton formed by Coulomb interaction appears 
and moves in the Cu02 plane. The two-peak structure 
of the observed spectra is attributable to the structure 
of the Cu3d partial density of states mixed with the 02p 
band, rather than some kind of bound exciton mode, (ii) 
The present formulation by the perturbative method is 
applicable to more complex and realistic electronic struc- 
tures, in contrast to small cluster calculations and exact 
diagonalization techniques, whose applicability is limited 
to simple solids consisting of a few kinds of elements. 
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